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Ph! Abstract 



In the generic supersymmetric standard model which had no global symmetry 
enforced by hand, lepton number violation is a natural consequence. Super- 
symmetry, hence, can be considered the source of experimentally demanded 
beyond standard model properties for the neutrinos. With an efficient for- 
mulation of the model, we perform a comprehensive detailed analysis of all 
(fermion-scalar) one-loop contributions to neutrino masses. 
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I. INTRODUCTION. 



Low-energy supersymmetry (SUSY) is the most popular candidate theory for physics 
beyond the Standard Model (SM). The most extensively studied version called the minimal 
super symmetric standard model (MSSM) has an extra ad hoc discrete symmetry, called R 
parity, imposed on the Lagrangian. It is defined in terms of baryon number, lepton number, 
and spin as, explicitly, R = (— i) 3B + L + 2S . The consequence is that the accidental symmetries 
of baryon number and lepton number in the SM are preserved, at the expense of making 
particles and superparticles having a categorically different quantum number, R parity. The 
latter is actually not the most effective discrete symmetry to control superparticle mediated 
proton decay [1], but is most restrictive in terms of what is admitted in the Lagrangian, or 
the superpotential alone. 

R parity also forbids neutrino masses in the supersymmetric SM. However, the recent 
data from the solar and atmospheric neutrino experiments can be interpreted in terms 
of massive neutrino oscillations. Thus, the strong experimental hints for the existence of 
(Majorana) neutrino masses [2] is an indication of lepton number violation, hence suggestive 
of R-parity violation. Giving up R parity, a tree-level neutrino mass can be generated 
through diagonalization of the neutrino-neutralino mass matrix. At the 1-loop level, all 
three neutrinos will become massive. There is then no need to introduce extra superfields 
beyond what is required by the SM itself to describe neutrino phenomenology. 

There is certainly no lack of studies on various RPV models in the literature. However, 
such models typically involve strong assumptions on the form of R-parity violation. In most 
cases, no clear statement on what motivates the assumptions taken is explicitly given. In fact, 
there are quite some confusing, or even plainly wrong, statements on the issues concerned. 
It is important to distinguish among the different RPV "theories", and, especially, between 
such a theory and the unique general supersymmetric standard model (GSSM) [3,4]. The 
latter is the complete theory of SUSY without R-parity, one which admits all the RPV terms 
without a priori bias. In the GSSM, RPV terms come in many different forms. In order 
not to miss any plausible RPV phenomenological features, it is important that all of the 
RPV parameters be taken into consideration. A clear listing and discussion of all these 
is recently presented in Refs. [5], under the framework of the single- VEV parametrization 
(SVP) [6,7]. The latter, summarized below, is an optimal choice of flavor bases that helps 
to guarantee a consistent and unambiguous treatment of all kind of admissible RPV terms 
with complete RPV effects on tree-level mass matrices for all states including scalars and 
fermions maintaining the simplest structure. Following the formulation, we present here a 
complete list of all the neutrino masses contributions up to 1-loop level. 

A (Majorana) neutrino mass term violates SM lepton number by two units. The ex- 
perimental evidence for neutrino masses comes in through indications of flavor oscillations, 
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which requires mass mixings of the flavor states, u e , and v T . Hence, we want neutrino 
mass terms that have lepton flavor violation (LFV). The latter is a generic consequence of 
R-parity violation. To put it in another way, the GSSM in fact contains many couplings 
that has one unit of LFV. Any combination of two of such couplings may be able to give 
rise to a neutrino mass term. Since the expected sub-eV neutrino masses are essentially the 
strongest source of upper bounds on such couplings [8] up to the present moment, we have 
no way to tell which particular combinations of couplings do saturate the bounds and give a 
dominant contribution to a neutrino mass term. In fact, each term also depends on a set of 
(R-parity conserving) SUSY or MSSM parameters. We do not have much knowledge on the 
SUSY parameters beyond some lower bounds on a set of related experimental parameters 
(mainly) from collider machines. In relation to the neutrino mass contributions discussed 
here, the set of SUSY parameters are typically taken as fixed by one generic SUSY mass 
scale. Changing the latter of course changes the actually neutrino masses resulted. More 
importantly, it is not totally clear whether some phenomenological hierarchy among values 
of the different SUSY parameters, may be together with some hierarchy among the values 
of the parameters with LFV, would not give a picture on the relative importance of the var- 
ious neutrino mass terms different from what one may expect from such the kind of highly 
simplified analyses. Thus, it would be useful to have a complete list of such neutrino mass 
terms without much a priori assumption involved. 

Guided by theoretical prejudices or otherwise, many different pieces of such neutrino 
mass terms have been studied [9,10] (see also Ref. [11] for a more updated list of references). 
More recently, there are attempts to give the more complete story. In particular, Ref. 
[12], gives the general formulae for neutrino mass contributions up to the full 1-loop level. 
However, the latter analysis is not formulated under the SVP and any detailed discussion is 
limited to a scenario where the "third generation couplings dominate" . Among the trilinear 
RPV couplings, this amounts to admitting only non-zero A'^'s and Aj 33 's, though all nonzero 
bilinear RPV are indeed included. The maximal mixing result from Super-Kamiokande may 
bring that wisdom of "third generation domination" under question. Refs. [10] and [13], 
for example, illustrate how no (family) hierarchy, or even an anti-hierarchy, among the 
RPV couplings may be preferred. More important to our perspective here, the study has 
assumption on the Bi parameters and is interested in the numerical study of a specific high- 
energy scenario. Here, we aim at a more detailed analytical study on the different pieces of 
contribution instead. With the help of a more simple theoretical framework, the SVP, we 
follow the basic approach of Ref. [12] and give a more transparent list of formulae, as well as 
pushing on to give much more detailed analytical results of each individual neutrino mass 
term. 

The basic approach of Ref. [12] is to give each 1-loop neutrino mass diagram in terms 
of effective couplings of the mass eigenstates of various scalars and fermions running inside 
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the loop, using a formula from the so-called "effective mixing matrix" method [14]. Details 
of all the admissible RPV contributions to all the scalar, as well as fermion, mass terms 
under the SVP framework are very manageable [16,4]. The complete expressions, together 
with useful perturbative diagonalization formulae for the interesting elements of the mixing 
matrices are listed in Ref. [4]. We use below exactly the same notation as presented in 
details in the latter reference, which is taken as the background of the present presentation. 
Our goal is to present the exact analytical expression for each neutrino mass term, and the 
approximate dominating result from each term under very mild assumption. The major 
part of the approximation is the perturbative diagonalization formulae of the mass matrices, 
which are well founded on the smallness of the neutrino masses. The approximation also 
helps to extract the major RPV parameter dependence of each mass term and, hence, is an 
important target of the present study. 

There is actually a detailed analysis of all the neutrino mass terms pretty much in the 
same spirit of present study published [11]. The latter reference also essentially adopted 
the SVP framework. However, mass insertion approximation is used to obtain the results 
based on the use of MSSM states. Our approach here may be a more direct and transparent 
alternative. Having results from both approaches also serving as a counter checking and 
helps to illustrate more clearly some of the subtle points involved. Ref. [11] also has some 
very different emphasis in their discussion. Hence, we consider the present study necessary to 
complete the story of neutrino masses in GSSM (or from R-parity violation). Moreover, our 
exact formulae, in terms of mass eigenstates running inside the loop allow direct numerical 
calculations of the neutrino mass results free from any approximation. We are also working 
on a detailed study of radiative neutrino decay within the model [19], to which the present 
paper also gives the necessary background. 

It should be emphasized here that it is not out intention to discuss scenarios within the 
general model that could fit the experimental date. There being such a large number of 
lepton number violating parameters within the GSSM, phenomenologically viable scenarios 
will not be difficult to find. The beauty of the GSSM in explaining the neutrino data 
is that the parameters responsible will also give a rich collections of other experimental 
signals. More studies of various aspects of the model, and constraints from various SUSY 
and LFV searches in the future may give much better guideline for picking the real interesting 
scenarios. The goal of the present study is to provide a useful better reference for such efforts. 

It should also be noted that we do not include here results of the gauge boson loop 
contributions. Such contributions have been studied (see for example article by Hempfling 
in Ref. [9]). They represent a small correction to the tree- level results which could be 
absoluted into a renormalization of the tree-level lepton flavor violating parameters, as also 
pointed out in Ref. [11]. Hence, we focus only on the finite 1-loop contributions that will give 
new structure to the neutrino mass matrix (same strategy was adopted in Ref. [11]). Namely, 
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we focus on the fermion-scalar loop. Our formulation in terms of mass eigenstates for the 
fermions and scalars inside the loop allow us to identify explicitly the Goldstone mode. Any 
Goldstone mode contribution is taken out from the summation over mass eigenstates running 
inside the loop. In fact, the Goldstone modes are of course unphysical and calculation of their 
contribution gauge dependent. The diagrams with the Goldstone modes only form gauge 
invariant sets with the corresponding gauge boson and ghost diagrams added together. If 
one wishes, one could take it that we are doing the fermion-scalar loop calculation in the 
unitary gauge. The gauge choice does not matter, as we are focusing on the sector of physical 
scalar bosons and fermions only. Including the Goldstone modes in our fermion-scalar loop 
calculations would rather be inconsistent. 

In Sec. II below, we give a brief summary of the basic formulation of GSSM used. Readers 
are referred to Ref. [4] for details. Sec. Ill then starts on the neutrino mass discussion. While 
the tree-level neutrino-neutralino mass matrix is quite well known, we present some of the 
details here for completeness. The presentation also sets the stage for the discussion of the 
1-loop contribution calculation. All the basics of the 1-loop analysis is presented in the 
latter parts of the section. The next section discusses some details of the results in the 
way outlined above. Some of the detailed listing of individual terms are, however, left to 
the Appendices. In Sec.V, we present a brief discussion on the application of the results to 
numerical studies, while any detailed numerical studies will be leave for future publications. 
Finally, we conclude the paper with some remarks in Sec. VI. 



II. BACKGROUND OF THE GSSM 



Let us start with summarizing our formulation and notations here; readers are referred to 
Ref. [4] for more details. The most general renormalizable superpotential with the spectrum 
of minimal superfields containing all the SM states can be written as 

W=e ab y a H a a L h a + KME h M + X ajk L a a Q)b c k + \ \ a , k L a a L h p E^ + \ X' l]k U^D^D c k , (1) 

where (a, b) are SU(2) indices, (i,j,k) are the usual family (flavor) indices (going from 
1 to 3). The 4 L Q 's, with the (a, (3) indices as extended flavor indices going from to 
3, include the usual leptonic doublets and the H d doublet. Four doublet superfields with 
the same quantum number are needed for gauge anomaly cancelation. The four are not 
a priori distinguishable. The rest of the superfield notations are obvious. Note that A 
is antisymmetric in the first two indices, as required by the SU(2) product rules, shown 
explicitly here with e 12 = — e 21 = 1. Similarly, A" is antisymmetric in the last two indices 
from SU(3) C , though color contents are not shown here. 

Doing phenomenological studies without specifying a choice of flavor bases is ambiguous. 
It is like doing SM quark physics with 18 complex Yukawa couplings, instead of the 10 real 
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physical parameters. As far as the SM itself is concerned, the extra 26 real parameters are 
simply redundant. There is simply no way to learn about the 36 real parameters of Yukawa 
couplings for the quarks in some generic flavor bases, so far as the SM is concerned. For 
instance, one can choose to write the SM quark Yukawa couplings such that the down-quark 
Yukawa couplings are diagonal, while the up-quark Yukawa coupling matrix is a product 
of (the conjugate of) the CKM and the diagonal quark masses, and the leptonic Yukawa 
couplings diagonal. Doing that has imposing no constraint or assumption onto the model. 
On the contrary, not fixing the flavor bases makes the connection between the parameters 
of the model and the phenomenological observables ambiguous. 

In the case of the GSSM, the choice of flavor basis among the 4 L^'s is a particularly 
subtle issue, because of the fact that they are superfields the scalar parts of which could 
bear VEVs. A parametrization called the single- VEV parametrization (SVP) has been 
advocated since Ref. [6]. The central idea is to pick a flavor basis such that only one 
among the L a 's, designated as Lq, bears a non-zero VEV. There is to say, the direction 
of the VEV, or the Higgs field Hj, is singled out in the four dimensional vector space 
spanned by the L^s. Explicitly, under the SVP, flavor bases are chosen such that : 1/ 
(Li) = 0, which implies L = H d ; 2/ y%(= X jk = -Xjok) = ^ diag{mj, m 2 , m 3 }; 3/ 
Vjk(= Xojk) = ^diag{m d ,m s ,m b }; 4/ yf k = ^\£ KM diag{m M , m c , m t }, where v = V2 (L ) 
and v u = \/2(H u ). A point to note is that the mi's above are, conceptually, not the 
charged lepton masses. The parametrization is optimal, apart from some minor redundancy 
in complex phases among the couplings. We simply assume all the admissible nonzero 
couplings within the SVP are generally complex. The big advantage of the SVP is that 
it gives the complete tree-level mass matrices of all the states (scalars and fermions) the 
simplest structure [4]. 

Following our notation above, the soft SUSY breaking terms of the Lagrangian, can be 
written as follow : 



V S o{t = £a b B a H^L b a + A 7 ^ Q1H b u Uj + A^H^Q b Dj + A E ij H d l L b Ef 



+ h.c. 



1 

2 

+ Q^rhl Q + &rf$ U + D^rhl D + DmfL + E^m 2 E E + ml \H,, 



A% k L1Q)D c k + - A\ jk LfL)E c k 



|2 



■ii | 



where we have used Hd in the place of the equivalent L among the trilinear A-terms. Note 
that L^mfL, unlike the other soft mass terms, is given by a 4 x 4 matrix. Comparing with 
the MSSM case, m. 2 corresponds to ml while m 2 's give new mass mixings. The other 

' L 00 1 H d L 0k ° ° 

notations are obvious. The writing of the soft terms in the above form makes identification 
of the scalar mass terms straight forward. Recall that only the doublets H u and bear 
VEVs. The A-terms in the second line of Eq.(2) hence do not contribute to scalar masses. 
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The SVP formulation also gives the complex equations 

Bi tan/5 = mj. + /x* , (3) 

reflecting the removed redundancy of parameters in a generic L a flavor basis. They are 
nothing but the vanishing tadpole equations. They give consistence conditions among the 
involved parameters that should not be overlooked. The equations suggest that the B^s 
are expected to be suppressed, with respect to the B , as the /Vs are, with respect to fj, . 
The m^. parameters in particular are missing in some of the relevant discussions in the 
literature. From a different perspective, one may tend to think that the parameters are 
similar to the rh%. parameters linked to soft flavor mixings. However, fixing m^. in Eq.(3) 
leads to definite relations between a Bi and a /Xj term, which may not be satisfied. The 
parameters Bi, fa, and m^. are not independent free parameters, because of the fact that 
freely chosen values of the set of parameters in a top-down approach, in general, do not land 
the model automatically into the single- VEV basis. The tadpole equations are incorporated 
completely into the scalar mass matrices involved in our calculations [4]. 

III. NEUTRINO MASSES 

The GSSM has seven neutral fermions corresponding to the three neutrinos and four, 
heavy, neutralinos. The heavy states are supposed to be mainly gauginos and higgsinos, but 
there is now admitted (RPV) mixings among all seven neutral electroweak states. In the case 
of small /ij's of interest, it is convenient to use an approximate seesaw block diagonalization 
to extract the effective neutrino mass matrix. Note that the effective neutrino mass here is 
actually written in a basis which is approximately the mass eigenstate basis of the charged 
leptons, i.e., the basis is roughly {v e , v^, u T ). The tree-level result is very well-known [9,10]. 

A. Getting the Neutrinos among the Neutral Fermions 

We use the basis (—iB, —iW, W u ,W d , Z" , Z° , Z° ) to write the 7x7 neutral fermion mass 
matrix Aiu- Note that W d = while fff is the charge conjugate of the higgsino W u . For small 
/ii's, we have ( ZJ, 1°, Z° ) ~ [y e , z/ M , v T ) [4]. The symmetric, but generally complex, matrix can 
be diagonalized by using unitary matrix X such that 

3CM„X = diag{M x o} . (4) 

Again, the first part of the mass eigenvalues, A^o for n = 1-4 here, gives the heavy states, 
i.e. neutralinos. The last part, il^o for n = 5-7, hence gives the physical neutrino masses. 
The mass matrix M. M can be written in the form of block submatrices: 
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M 



m,. 



(5) 



where M. n is the upper- left 4x4 neutralino mass matrix, £ is the 3x4 block, and m° is the 
lower-right 3x3 neutrino block in the 7x7 matrix. In the interest of small neutrino masses, 
a perturbative (seesaw) block diagonalization can be applied. Explicitly, the diagonalizing 
matrix can be written approximately as 

The tree-level effective neutrino mass matrix may then be obtained as 

-cm; 1 er M n (ac f ) = - e ac f 

M| cos 2 /3 (M cos 2 6^ + M 2 sin 2 6U) 



(6) 



det(A(, 



( IH ) 



where 



det(Af„) = fio \-iJi M 1 M 2 + Ml sin2/3 (M 1 cos 2 6U + M 2 sin 2 ^) 

is equivalent in expression to the determinant of the MSSM neutralino mass matrix. 
It is obvious that the 3x3 matrix ( fa fj,j ) has only one nonzero eigenvalue given by 



(7) 



(8) 



fi B = \fa\ + \fi 2 \ + \fi a 



(9) 



We can define 



Rr — 



El 

M5 

El 

Ms 

El 

M5 





_M3_ 



Vlwl 2 +lMl 2 



M5 

mi Ei 



M5 Vl^| 2 + lM3| 2 

Ml Ms 



Vl^l 2 +lwl 2 msVI^P+ImsI 2 / 
Then, we have iZjT ( /i« /ij ) i? 5 = diag{ ^ 2 , 0, }. Here, /i 5 and 



/'-2 



A*3 



(10) 



are taken as real 



and positive. With this result, we can write the overall diagonalizing matrix X in the form 



X ~ 



-(ACf) 





-(ACf)^ 



e* (Al; 1 f ) i? 5 



(11) 



where i?„ is a 4 x 4 matrix with elements all expected to be of order 1, basically the 
diagonalizing matrix for the Af n block and e lC - is a constant phase factor put in to absorb 
the overall phase in the constant factor in the expression of Eq.(7) so that the resulted 
neutrino mass eigenvalue would be real and positive. The matrix X contains the important 
information of the gaugino and Higgsino contents of the physical neutrinos. This is given 
by the mixing elements in the off-diagonal blocks. The Z matrix in itself gives similar 
information for the effective SM neutrinos (flavor states). The latter matrix may be more 
useful in the analysis of neutrino phenomenology. 
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B. Approach to 1-loop Neutrino Masses Calculations 



Following Ref. [12], we use the 1-loop (renormalized) mass formula from the "effective 
mixing matrix" approach, giving a fermion mass matrix as 



M/V) =M(Q) + rV) ~-[m(q) s( p 2 ) 



S(p 2 )-M(Q) 



(12) 



Note that M(Q) is the DR renormalized tree- level mass (matrix), while II and S the 
contributions from one-loop self-energy diagrams with and without chirality flip. We have 



AVV) 



r 





5£ 5(m° u ) 



(P 2 ) , 



(13) 



where 



SM n ( P 2 ) 


= n„( P 2 ) - 


1 r 

2 - 


sap 2 ) 


= n,(p 2 ) - 


1 r 

2 - 


s<(p 2 ) 


= ixV) - 


1 r 

2 - 



M n T, n {p 2 ) + ^n{p 2 )M, 



(14) 

with the explicit renormalization scale (Q-) dependence of the tree-level parameters dropped. 
Seesaw diagonalization of Mm {1) yields the 1-loop result 

W (1) - -e -m; 1 f + *«) - ^ aC f - e .m; 1 *f + e a*; 1 sm. a-c f 



+ i e^ac f + ^ aC f K + e >C n n ac f , 



(15) 



where we have dropped the p 2 dependence. As discussed below, the p 2 should be taken as 
at the scale of the mass M. u itself. Hence, in the application here to calculate the neutrino 
masses, the p 2 in (m u ) w (p 2 ) may be taken as practically zero. An important point to note 
here is that the and S n terms all cancel out and disappear from our final result for (m„) (1) . 
We refer the reader to Refs. [12,14] for further discussion on the merits of the approach and 
references to related works. 

At this point, some remarks on the renormalization issue are in order. The issue has 
been well-addressed in the papers by Hempfling and Hirsch et.al. [9] on calculations starting 
from tree-level mass eigenstates. The analog formula to Eq.(12) above is, under the DR 
scheme, 



M polc = M DR {Q) + AM(p, Q) 
where the 1-loop correction part is given as 



(16) 
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AM(p,Q) = [AM(p)] A=0 , 



(17) 



i.e. the two-point functions involved are calculated by subtracting the term proportional 
to the regulator A = — j E + In Air of dimensional reduction. There is some ambiguity 
in the choice of Q in the evaluation of the off-diagonal two-point functions. As pointed 
out in Hempfling's paper, the effect of the ambiguity is of higher order. In the "effective 
mixing matrix" approach [14], the equation is casted in the electroweak state basis instead 
to arrive at Eq.(12), which upon seesaw block-diagonalization yields Eq.(15). Now, p 2 is 
practically zero, as we are calculating only diagrams with neutrinos on the external legs of 
the two-point functions. The rest are only tree-level mass matrix entries (to M^) coming 
into the formula as mixing matirx elements between the neutralino and neutrino blocks. 
The result for (m„) (1) is however Q-dependent. Apart from the Q-dependence in the former 
set of parameters, there is also the Q-dependence coming up from the calculation of the 
two-point functions under the DR scheme as in Eq.(17). Furthermore, there are the full set 
of couplings involved in such calculations, which should be taken as running couplings at 
the scale Q. In the straight formal sense, the pole mass formula gives result M polc that is 
Q-independent. However, in the application to obtain M. u m {p 2 ) and its subsequent use in 
any explicit calculations, some residual Q-dependence is difficult to avoid. 

Since we are interested in radiative neutrino mass generation from superparticles, we 
may take Q as M SUSY , or roughly the electroweak scale. Below the scale, the superparticles 
decouple and the neutrino mass terms can only be expressed in terms of five-dimensional 
operators of the SM. Strictly speaking, one get the correct pole mass for the neutrinos only by 
running the operators to the neutrino mass scale through the corresponding renormalization 
group equations. However, such effects are minimal [15]. Apart from yielding the neutrino 
mass matrix in the more interesting flavor basis, Eq.(15) also avoids the superficial singularity 
reflecting the arbitrariness in the diagonalization of the mass matrix with degenerate massless 
neutrinos at tree-level. Furthering, the MNS (neutrino mixing) matrix obtained from the 
diagonalization of (m u ) (1> maintains a full unitary matrix. 

The full 7x7 neutral fermion mass matrix has four heavy and three very light mass 
eigenvalues, corresponding to the neutralinos and neutrinos; and we are essentially only 
interested in the neutrino states. For a general mass calculation at 1-loop, we must choose a 
renormalization prescription for each of the tree-level parameters appearing in the full mass 
matrix M. U) and will have to worry about the renormalization scale dependence issues of 
such parameters. The "effective mixing matrix" formula [cf. Eq.(15)] avoid the complication 
as the neutrino mass results depend only on p 2 explicitly, which is practical zero, as pointed 
out in Ref. [12]. The parameters involved in the neutrino mass generation are then taken as 
running parameter at the scale of interest. Renormalization scale dependence comes in only 
through the E„ part (readers of refer to Ref. [14] for more details) and the effect is small. 
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The T, u part itself is mostly not very important, as can be easily seen in Eq.(15). 

Our neutrino mass formula [Eq.(15)] calls for a seesaw type block diagonalization of 
the mass matrix M. M up to 1-loop order. The diagonalizing transformation corresponds to 
the matrix Z of expression (6). The tree-level contribution, given by the first term in the 
formula, is obviously seesaw suppressed (by the neutralino mass scale). The second term 

gives the direct 1-loop contributions. However, there are parts of U u that involved other 
suppression beyond the loop factor. A typical example is the pure gaugino loop, or GH-loop 
[20], diagram contribution which can be interpreted as requiring seesaw induced Majorana- 
like "sneutrino" mass to give a nonvanishing result [5]. They may be called pseudo-direct 
1-loop contributions. For the rest of the terms in Eq.(15), are indirect 1-loop contributions, 
which has part of the basic seesaw suppression going along. These include results from 1- 
loop diagrams contributing to the off-diagonal blocks of the M. M matrix, from H v diagrams, 
as well as from diagrams contributing to the diagonal block Ai n . The last one, given by 
the last term in the formula, gives no interesting features. It can be absorbed, for instance, 
into the tree-level result (first term) by replacing Ai n there with the 1-loop corrected result. 
And, from the related calculations within MSSM, we know that the correction is about 
6% [14]. In fact, the flavor conserving part of the contributions involving T, v is similarly 
uninteresting. However, the part of the latter with LFV may be of interest. 

To calculate explicitly the various neutrino mass contributions using the above formula, 
we need to have the effective couplings of the electroweak state neutral fermions to possible 
scalar and fermion mass eigenstates running in the quantum loop. The neutral fermion 
themselves, together with the nine neutral scalars of the model, give a class of neutral loop 
contributions. Obviously, the loop with the neutralino states dominates here. The effective 
couplings, to be given below, involve diagonalizing matrix elements of the states contributing 
to the states running inside the loop. For the fermion part, it is the X matrix discussed 
above. Similar perturbative diagonalization expressions for all the other matrices, those for 
the charged fermion, charged scalar, down-squarks, as well as the neutral scalar sector are 
discussed in details in Ref . [4] . We refrain from repeating the long list of such formulae in this 
paper. Most parts of the notation used, as will appear below, are quite easy to appreciate. 
Readers interested in checking any details on the derivations of the results, however, would 
need to use Ref. [4] extensively. 



For the neutral loop contributions, we start with the effective interaction for the external 
Z°'s with internal mass eigenstates, 



C. Neutral Loop Contributions 




2 



1+75 

2 



*(x°) C + h-c. 



(18) 
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where |(1 =F 7s) are the L- and f?-handed projections. We have 

A(l = \ [tan^X,; - J£,*J [P ( s l+2)m + ^ £f i+7) J , (19) 
and A/" L = Af R * . The direct 1-loop contributions is given by 1 

IT = AP* A/"** Af o i3o(V , M 2 , M 2 ) , (20) 

j 87r sin 2 # mm ]nm ^ x ™ ' 

where the loop function Bo is defined in the limit of p 2 — > by 

Sob 2 , < ml) = -^-r ^ 4 - ^ ^ + 1 (21) 

As will be shown explicitly below, this result is the gauge loop contribution first discussed 
in Ref. [20]. Note that X is the matrix that diagonalizes the seven neutral fermions, as 
discussed explicitly above. Among the seven fermion (tree-level) mass eigenstates denoted 
by the sum over n here, contributions from the n = 5-7 states are certainly negligible. The 
sum of m runs through the nine physical neutral scalar states. The states, together with 
the unphysical Goldstone mode, are obtained from the 10 x 10 neutral scalar mass-squared 
matrix to be diagonalized by V s . We refer readers to Ref. [4] for details on the scalar sector. 
The set of coupling vertices may also be combined to give contributions to the self energy 
function E„. We have 

K = HU^L ®{v\ Ml, Ml) , (22) 

where the loop function B\ is defined by in the limit of p 2 — > by 

1 



Bi(p 2 ,m\,m 2 2 ) = | 




2 j m\ 1 / m\ + m 2 . " 
m\ — m" 2 .) m 2 . 2 I m 2 — m 2 . 



(23) 



Here, we have all fermions involved being Majorana fermions. We compose the 4-spinor \I> by 
* = ( ^ I wnere we nave Vfe = —i^^L- A mass term has = ^ip R + 1$ tp L . For instance, 

the if$ ip L part can be written as ^ ia 2 ip L . The A£^, vertex brings in a ip L . while a matching A/^ m 
vertex brings in a tp R . = ia 2 . With the proper handling of the fermion wavefunction, 11^. has 
contributions proportional to A(^ m A^ m which is equivalent to the M^ m M^ m used here. Just for 
the 1° states, the use of A(f m can only be seen as the intrinsic left-handed nature of the states. 
However, the explicit use of their charge conjugate to compose the 4-spinors and derive the effective 
couplings is necessary to complete the formulation, say in the case of the C£ m for the charged loop 
discussed below. 
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For the indirect 1-loop contributions, we need 

KL = \ [tan« - X? n ] [V* m + i V°J , (24) 

ML = -\ [tan^; - X£\ [Vl m + i V° m ] , (25) 

M£nm = 2 ^3*n [Am + * ^6m] ~~ 2 [^ > (2+a)m ~~ * ^(7+a)m] > (26) 

M£ m = -\ tan^X 3 ; [Z£ m + i P 6 S J + 1 tan^X ( : +Q)n [Df 2+a)m - i Vf 7+a)m ] . (27) 

The list of extra M R * terms each combines with the M R * m to give a neutral loop contribution 
to IT ? . 



D. Charged Loop Contributions 

The effective interaction for the external 1° with (colorless) charged fermions and scalars 
inside the loop is given by 



r L - + C 1 



^(Xn)^ +h-C. , (28) 



where 



7/ A* 
/ifl _ tki 17 n'* 2%kh\r T)V 
— Vfi+2W M>™ Vf/,.+oW M 



mm 



(i+2)n ^2m — y (/i+2)n ^(fc+2)m 



Cinm - -Uln ^(i+2)m + ^ ^2n ^ ) (i+5)m ~ £fyi+2)n ^(fc+5)m ■ ( 29 ) 

Here, V^A^c — diagiM -} where is the 5x5 charged fermion mass matrix. Matrix 
D l diagonalizes the mass-squared matrix of eight scalars of unit negative charge (see Ref. [4] 
for details). The latter includes again the unphysical Goldstone mode to be dropped from 
the sum over m. 

Charged fermion loop contribution to direct 1-loop neutrino mass could then be easily 
obtained as 

= C r* c l M _ B ( 2 M 2 M 2\ , ii — • ) (30) 

ij 8ir sin 2 6U mm Jnm Xn ^ x ™ ' 



L 

jnm. 



Unlike the case for the neutral loop result, the IT^ matrix written through the C H * C: 
co up led- vertices is not symmetric with respect to % and j. Hence, an explicit symmetriza- 
tion has to be performed, as indicated above. The symmetrization also takes care of the 
asymmetry with respect to L and R, automatically. Similarly, for the T, u part, we have 

Yf v = — \C L C L * +C R *C R } tf/p 2 , M 2 _ , M 2 ) (i< — >j). (31) 
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To go on to discussions of the indirect 1-loop contributions, we need the corresponding 
expressions of the C££ for the other four neutral fermions. These are given as follows, with 
obvious notations, 

C nm = -Uln Am - -f- U( k+2 )n ^(k+5)m ! (32) 
Ham — y ln ^lm , 

<*L = ; (33) 

-1 

V2 

l 1 r * * 1 

Cfrnm = ^ [U2n'D 2jn + U( k+ 2)n'D(k+2)m\ 5 (34) 

— t&n9 w 



Cwnm ~ /pr ^2n Am > 



C R 

Bnm \/2 

L tan^v 

Bnm \f^2i 



V2n T^lm + 2 V(k +2 )n ^ > (jfc+5)m 



^2n I*2m + ^(fc+2)n ^(fc+2)m I ■ (35) 



Combining a C fl * with a C L gives half of the charged fermion loop contribution, to the 
corresponding mass term; the other half is given by flipping L and R. For instance, the 
l°-B mass term, or 11^, is given by substituting C^ m for C L nm in Eq.(30), i.e., by a C R * m C~ him 
combination, as well as the combination C£* C L . 

' Brim trim 

There is also another type of contributions, namely the quark-squark loops. The direct 
1-loop part of such contributions is among the most well discussed. We summarize them 
here, under our notation, for completeness. We have 

K„ = ~ ^f: a C: C m dn B Q ( P \ ml Ml) ( i < > j ) (36) 



8vr sin 2 6U 



where 



nIR 


A'* 

ikn 
92 


U km > 


nlL 


A' k 

92 


T) d * 

U (k+3)m 



(37) 

and V d diagonalizes the 6x6 squark mass-squared matrix The structure is to be 

compared directly with those from the A-couplings above. For E„, we have 

^ = — \C L C IL * + C IR 'C IR } B X U,M} ,M?) (i< — > j) . (38) 
For the indirect 1-loop part, we need 
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nIR Vdn <T-\ 

Hhm „ km ' 



Onm — U (k + 3)m 1 l dy J 

9i 



C nm = ^L; (40) 

C£ m = -^t a n9 w Vf m . (41) 

We get the indirect 1-loop contributions by combining C' R(L> * with C'^^p, in the same way 
as we do in the above case of (colorless) charged fermion loop. 



IV. MORE DETAILED ANALYTICAL RESULTS 



In this section, we give more explicit details of the neutrino mass terms obtained by 
applying the formulae in the previous section. We list the result from different combinations 
of interaction vertices and go on to illustrate the content of these exact mass eigenstate results 
by extracting the dominating piece (s) in the mass eigenstate double sum. There, we give the 
"approximate" analytical results through the use of perturbative diagonalization expressions 
[4] for the elements of the various mixing matrices. Such perturbative diagonalizations have 
been illustrated to be very good approximations, which also serve to illustrate well the role 
of the various lepton flavor violating (LFV) couplings involved (see Refs. [17,18] for other 
illustrations). 

Note that we focus our discussions below only on the parts of the results that 
are particularly interesting to our analytical study. For instance, in the 11^ = 
- a a ™ 2a N R * M R * Mo B (p 2 , M\ , M„ 2 ) term, we focus on the N R * M R * Mo part. That 
is, we will drop the common pre-factor 8?r "^2^ and the loop integral £>o from all the neutrino 
mass term results given below. The following discussions do not include the part. The 
results of the latter are left all to an appendix at the end. They are included here mainly 
for completeness. It does not look like there is any important off-diagonal contribution, 
while diagonal contributions, as discussed above, only represent a universal correction to 
the tree- level result. 

A. Results for n£ and (n^M^f)^ 

• The result here may be written in the form of a single term as 
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- [tan^X ln - X 2n ] 2 A^o [D* i+2 -) m - i2?( S i+7 )J [^ > (j+2)m ~~ ^o+7)m] 



SiSj tan 2 /? 



- (tan^X ln — X 2n ) 2 



(n=l-4 dominate) . (42) 



The scalar sum is dominated by m — 1,2 and 7 contributions. We illustrate here 
only the dependence on the B { parameters and tan/3, with M s denoting a generic mass 
parameter at the slepton scale. 2 Note that we write the final result in the form such 
that the square-bracket [ ] contains a factor of order 1 (a pattern we stick to below), 
so that the reader can have an idea on the major parameters (those before the square- 
bracket) affecting the scale of the neutrino mass. The resultant proportionality of 
the mass term here to the product B; L Bj has been addressed and interpreted as the 
necessity for a Majorana-like scalar mass insertion to complete the diagram, in terms 
of complex scalars. When one follows such an interpretation to consider the scalar 
inside the loop as complex field with mass insertions put in on the line explicitly (as 
shown in Fig. 6 of Ref. [5] for example), a proportionality on B^Bj would likewise be 
resulted. The different m pieces in the scalar sum, however, cannot be put together 
at this level. Each piece involves actually a different value for M s and a different loop 
integral from a physical scalar of different mass running in the loop. In fact, if one 
naively takes a sum over m without considering the loop integrals, a zero result would 
be obtained for any T>^ m T>^ m with a ^ b. The lack of degeneracy among the scalar 
mass eigenstates is what makes a nonzero result possible. This is a common feature 
for the type of diagrams (see also Ref. [18]). Interestingly enough, for the present 
case under discussion, a pairwise degeneracy among the "scalar" and "pseudoscalar" 
parts of a complex scalar is enough to guarantee a null result. This is equivalent to 
the statement that the neutrino mass contribution is proportional to a Majorana-like 
mass term. It is illustrated here in our expressions as a consequence of the cancelation 
between V( %+2)m D( J+2)m and i 2 V( i+7)m D( j+7)m as well as between V( i+2)m D^ +7)m and 



2 Note that we have [£>^ +2 )i — ^(i+7)il — ~^m^ ~^ [tan/3 sina — cosa] = — [tan/3 sina — 

cosa] from our perturbative formulae on the V s elements. One may also check the other pieces. 
Take the m = 2 piece, for example, we have then [^ +2 )i ~~ ^(l+7)il — ~~ W 1 [^ an fl cosa + sina]; 
for m = 7, [£f i+2 )7 - iV (i +7 ) 7 } - ^co^A^^ = i % ^3- The extra factor of i guarantees a 
cancelation with the m = 1 and 2 terms if the (m = 7) "pseudoscalar" is mass degenerate with 
the latter "scalars", as [tan/3 sina — cosa] 2 + [tan/3 cosa + sina] 2 = co ^^ . Hence, to illustrate 
the generic result, we write the dominating [^ i+2 ^ m — *^(i+7) m ] result as ~' B ^f n/3 . This is used 
throughout the section. 
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/ D(i+7)m'D(j+2)m f° r eacn single m value, from our perturbative expressions for the mixing 
matrix elements. 

Next, we come to the (lit .M" 1 part. The dominating results from all the individual 
terms of the form have a common proportionality to the combination of LFV parameters 

Bi fij (tan/3) . 

Again, the contribution mainly involves diagrams with a (physical) neutralino, together with 
a neutral scalar, running in the loop. As noted above, the lack of mass degeneracy among 
the scalars is essential for a nontrivial result. Note that upon the necessary symmetrization 
not explicitly shown, we will have also the Bj fa (tan/3) parameter combination coming in. 

All the different terms in this class have very similar structure. We discuss here only an 
illustrative term, and leave the rest to Appendix A below. Let us take a look at the term 
^(.M; 1 ^, It is given as 

[i„M 1 M 2 - M z 2 sin/3 cos/3 (M 1 cos 2 6w + M 2 sm 2 ^) 1 2 
-Hi i Nn^X ln - X 2n ] M# 

' \P{i+2)m ~ i £fi+7)nJ Km - * ^Tml 



Bi 11 j tan/3 



(ioM,Jli-M z 2 sin/3 cos/3 (M cos 2 6» w + M 2 sin 2 (9 w ) 1 2 



(43) 



Note that from the general flavor structure of the model, one expect IT^ to have an expres- 
sion similar to H"^. above with index j replaced by a 0, i.e., IT^ ~ tan 2 /3 | (tan^Xi„ — 
X 2n ) 2 A^o . Observing that is a parameter of the same order as the generic mass scale 
parameter il£, 2 tan/3 we do see an agreement here. 



B. Results for the Charged and Color Loops 



Similar to the neutral loop case above, each term in the charged loop contributions to 
the IPs has a scalar part involving D l am T>l m which would give a null result for a ^ b if 
summed over m naively. The different loop integrals from the lack of scalar degeneracy is 
what guarantees nontrivial results. The fermionic part is more interesting. For illustrative 
purpose, we take an expression of the form V* i+2 ) n M- U\ n . Here, n = 1 and 2 give the 
chargino state contributions, with a large M- mass but a more suppressed V(* +2 ) n mixing. 
The results are given by -jgy m; /i« U\ n and -jf^- rrii fa Um, respectively. On the other hand, 
the n = i + 2 term involves a small fermion mass M - — but a less suppressed mixing 

x i + 2 

of t/i(j + 2) — ilJ^w^i ^ Dropping all the factors of order 1, we have all three terms giving 
contribution of roughly the same order of magnitude, all proportional to ^jf 1 , where M s 
again denotes a SUSY scale mass parameter here corresponds, more exactly, to a chargino 
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mass. This kind of feature is quite common in the charged loop results below. We illustrate 
results by dropping all the order 1 parameters and using the generic mass parameter M s 
representing chargino as well as slepton mass scale. 

There are six terms to the U.^.. result. We mostly just list them, while drawing attention 
to particularly interesting features. Note that the necessary symmetrization is not shown 
explicitly. 

Vh v * m tj n i v i* „ Vh mjHBj tan/3 
" — » / ( i +2)n U ln V 2m D (j+2)rn ^— — 3 (44) 

i/2 i/2 1V1 S 



The scalar part result here is mainly from — 

i/2 i/2 i/2 i/2 1V1 S 



(45) 



fij tan/3 



The scalar part result here is mainly from -D 2 Z (5+j) c: 

^ V*= v * M t t n i W* „ Vh Kk m m k (i k tan/3 

" — — " V(i+2)n U (h+2)n Am Z>(k+5)m ~ ~ ~ ~ ( 46 ) 

i/2 i/2 i/2 i/2 1V1 S 

Here, the result is from n = i + 2 which is interesting only ai h = i; hence, only that 
is shown in the sum over h. It is the SUSY analog of the Zee diagram, discussed in 
Refs. [13,5]. For h ^ i parts, the result is much further suppressed (by another 
factor). The scalar part result is the same as the previous case. 

^ V ( * h+2)n M- U ln D{ k+2)m V[ ]+2)m - ^ ^ (symmetrization !) (47) 

We note here that the result is actually very sensitive to the i <-> j symmetrization. 
The dominant result in the expression above is from the case with the (j+2)th charged 
scalar running in the loop. This is approximately the l~ slepton. The symmetrization 
and the fact that = —Xjih suggest a perfect cancelation of the result in the limit 
of degenerate sleptons which correspond roughly to the l~ and l~ states. This has also 
been discussed in some detail in Ref. [5]. 

Ikj \kh - « M rr n l T> 1 * ~ ^ Xikh m ^ ^ ( M L)*jk , 

~J 2 ~^ V ^ n M - U2n D ^ m V ^ m ~~J 2 ^ 2 Ml (48) 

where 

i Ta2 \ r a* m r \/2 M w cos f3 „* \pl M w sin/3 

(M 2 RL ) jk = [A* e - fji tan/9] m 3 6 kj + 5A% ( ^ ) . (49) 

92 92 

gives the complete LR mixing of lj and l k states. The last part of the latter is a 
contribution beyond the well known MSSM parts. 
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Kkh Ajgrp J Ajfcfe Xjhp m h {M\ L )* k 

— — ^(h+2)n M x; ^(9+2)" D (fc+2) m %+5) m - — ~ 77^ (50) 

m i/2 i/2 i/2 J «j 

This is the most well known part of II£ result discussed extensively in the literature. 
Note again the extra (last) term in the LR mixing (M.^ L )* k . Its contribution to neu- 
trino masses in the case of p ^ k may be particularly interesting. 



For the (II^At^ 1 ^ T )£- part, we present the long list of terms in Appendix B. In the 
neutral loop counterpart above, we see that the class of indirect 1-loop result all involve 
the combination Bi fij tan/3. Here for the charged loop results, we see the same parameter 
combination does give some important terms, but without the tan/3 factor. These are labeled 
as IT^o (M~n ^)3j — part 1 and n^ (M'^ — part 5 [with corresponding Il^ (M'^ ^)2j 
part] inside the appendix. In factor, these terms could easily dominate over the direct 1-loop 
terms from U.^ over. They provide neutrino mass contributions of order -^r- 

Another type of interesting terms are given by those labeled as (A / (~ 1 ^ T )ij — part 
4 [again with corresponding IL° fi {M.' r } £ r )2j part] and (.M" 1 — part 10 inside the 
appendix. We have, roughly, the results ^ or ^ 

Most of the other terms are actually not very interesting. They typically involve further 
suppression from factors such as ff-^-- However, one should note that for the large tan/3 
case, the i = 3 part has an order 1 coupling (essentially the r Yukawa) which renders the 
suppression not too strong. A careful numerical study will be necessary to check if there 
could be a scenario where such term could play a role. 

The quark-squark loop results are much more simple as a class. In fact, parallel structure 
between the \'^ k coupling terms and the A^ coupling terms can also be used to write down 
the results directly. In particular, for the indirect 1-loop part, we expected ^ 3 ^2 h or 
^ m^ H to ma tch the similar terms just discussed above. We list the details in Appendix C. 



V. CONCLUDING REMARKS 



We have listed and discussed the detailed results of all the neutrino mass terms within 
the GSSM, up to 1-loop order. Our approach gives expression for exact results, each to 
be obtained through a double summation over the fermion and scalar mass eigenstates 
running inside the loop. We further give approximate expressions of each of these terms 
through extracting the dominating pieces within the double summation and approximating 
the elements of the mass mixing matrices by perturbative diagonalization formulae. The 
validity of such perturbative diagonalizations are well founded on the experimental smallness 
of effects involving lepton flavor violation or R-parity violation. However, there are partial 
cancellations among pieces within the sum — - a result of a GIM type unitary cancellation 
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also pointed out in Refs. [17,18], rendering the approximate formulae agree only at order 
of magnitude level with the summed exact results. The latter is also cross-checked through 
numerical calculations, part of which is given in Appendix E for illustrative purposes. We 
most probably have given the results in more details than necessary. However, we emphasize 
that our ignorance about the nature of SUSY parameters, R-parity violating or otherwise, 
says that imposing much theoretical prejudice on the likely importance on some contribution 
over the others may be unwise. The detailed listing here is intended to provide a reference 
to later studies on any plausibly interesting scenario out the model. Numerical studies of 
the latter will be published independently. 
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APPENDIX A: DETAILS OF (L^A^f )~ TERMS 



n^(ACf)4; 



HoM^Mz - M z 2 sin/3 cos/3 (M 1 cos 2 w + M 2 sin 2 w ) 1 



det(M n ) 



[ta.n0 w X ln - X 2n ] 2 M x o 



Bi /ij tan/3 
M 2 



ftMMj - M| sin/3 cos/3 (M 1 cos 2 8 w + M 2 sin 2 6W) 1 
dct(M n ) 4 



(tan0 w X ln -X 2n ) 2 M x o 



■ (Al) 



This is exactly expression (43) which we repeat. 



M 2 cos 2 /3 (M cos 2 ^ + M 2 sin 2 w ) 1 2 
H deT^MJ 4 t ta n^X ln - X 2n ] M x , 



[£fi+2)m ~ i £fi+7)m] [Am _ 1 % 



-Bi fij tan/3 
Ml 



M| cos 2 /3 (M x cos 2 w + M 2 sin 2 w ) 1 



det{M T . 



- (tan^X ln - X 2n ) 2 M x o 



(A2) 
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Again, the H"^ term has a structure similar to that of n^. (or n^) with the replace- 
ment ofh d (= l ) by fy. 



n^WfH-- Parti 

M z cos/3 fi M 1 cos9 w 1 



det(M„) 



- X, n [tan^X ln - X 2n ] M x c [P ( s l+2)m - i £f i+7) J [D{ m + i V° m ] 



-Bi Hj tan/3 



M z cos/3 /UoM cos6* w 1 
det(M n ) 4 



X„ (tan0 w Xi„ - X 2 „) M x o 



(A3) 



(M; 1 f )2j part 2 

M z cos/3 /j, M 1 cosQw 1 



det(7W„) 



[ ;Z3 (j+2)m * I) fj+7)m] [^(2+Q)m + 1 ^(7+a)m] 



/Uj tan/3 



M z cos/3 /XoMi cos0 w 1 
&et(M n ) 4 



Xi„ (tan^X^ - X 2n ) M x o 



(A4) 



Here, we have different terms for a = 0-3, among which we show only the a = result. 
The a = 1-3 cases have obvious extra suppressions from the X^ +a ) n matrix element 
and a smaller scalar mixing part. The former has an extra -j^- factor while the latter 

introduces a — ^ factor. The overall H"^ results are not too different from the 

previous ones above either. 



^(ACOii-partl 

M z cos/3 ^ M 2 sin^, 1 



det(M„) 



Si /Uj tan/3 



(A5) 



M z cos/3 /x M 2 sin9 w 1 

det ^ ^ j tan# w J*3„ (tan0 w Ai„ - A 2 „) M x o 



n^(A4; 1 f) lj -part2 

M z cos/3 /j, M 2 sm6 w 1 



det(M r 



- t&n9 w X(4 +a ) n [t&n8 w X ln - X 2n ] M,o 



Bi Hj tan/3 



M z cos/3 /z M 2 sin6W 1 

dct ^ j ^ tan6» w JQ„ (tan^,X ln - X 2 „) M x o 



(A6) 



If one rotates the bino and wino into a photino and a zino, the photino would of course 
be decoupled from mass mixings with the neutral fermions. The closely related struc- 
tures of 11^ (M'n ^ r )2j and U"g {M~n reflect on that. One can certainly write 
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the two part of the results together through a U". 2 term with diagrams involving the 
zino part only. However, to the extent that photino and zino are not mass eigenstates, 
there is really not much to gain. 



APPENDIX B: DETAILS OF (U^Mr^^) TERMS 

• Kp (AC > 

Here, we introduce the order 1 constant 

yUoMMj - Ml sin/? cos/3 (M 1 cos 2 9 w + M 2 sin 2 6W) . , /T31 N 

°4 = TT771 \ Ms l m J 

det(M n ) 

to simplify the expressions, given as follow. 
Kp (AC e )4j Part 1 

^ ^ J 2 V ^ n *» C/ln Am ° 2m - J 2 ~Mf~ 64 (B2) 



n^ (AC f )<j - part 2 

j ~M s 77 (l+2)n x ' n {k+2)n 2m V ( k + 5 ) m \yj M] 4 ^ } 

Kp (AC eh part 3 

C 4 A ifch , A lfch //j m fe g| tan/3 

_ ^ ^ ^ V+2)n M x n U ^ D {k+2)m ^m--- ^ C 4 (B4) 

n^ (AC f ) 4j - part 4 

~* W s ~^J 2 V (W)nK- U ^)n D (k+2)m D (p+5)m _ -- — — — C 4 (B5) 

n^ (AC f )<* - Part 5 

Jf g — V (k+2)n M x - U ln V (k+2)m V (i+2)m M 2 °4 (BbJ 

n^ (AC f )<i - Part 6 

(M 2 RL y ik n 

N M.lh~ih {k+2)n x; (fc+2)m (J+5)m " 7 7 Mj 64 (B7) 
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n& {Ml 1 f) 4j Part 7 

~ N ~m~ s ~lhlh (p+2)n * ; (,l+2)n (p +2 > m < fc+5 ) ra #7 m| 4 ^ j 

Here, we need to use, in addition to above, expressions for the elements of the mixing 
matrix D[^ i+2 ^ ~ and D[^ i+ ^ ~ and introduce the order 1 constant 

Mj cos 2 /3 (M cos 2 ^ + M 2 sin 2 6W) , 
det(.M n ) 

to simplify the expressions. We also use M c to denote a mass parameter of the (phys- 
ical) chargino mass scale. The results are as follow 

Kh W £%" - P art 1 

f V - M x; **» D L DH +2)m * (BIO) 



m; ^ ^ M - ^ A - " ~^ n^r 3 ^ ( } 

n ^ (^n - P ar t 3 

H W s ln *» (h+2)n lm ^( fe + 5 )™ 3 ( ^ 

(- M « and (M; 1 f)y > 

We have noted above in the case of IT^ the close similarity between {M-'n £ r )2j 
and n^ u ^ (-M" 1 The story in the same here, between (Ai'^ e) 2 j and 

^v-b {M~n ^ T )ij) with some exception. Note that from comparing Eqs.(34) and (35), 
we can see that there is and extra term in C^ nm without the matching partner in C£ nm . 
We list below all results of the bino case, namely, {M-'n Oir Among the ten 
parts listed below, 1 to 7 have the wino counterparts in {M~^ £?)2j, to be given 
with an extra factor ta ^l j^, which we do not list explicitly. Parts 8 to 10 have no 
wino counterparts. We also introduce order 1 constants 

= tanfly M z cos/3 fx M 2 sin8 w 
l ~ ^2 det{M n ) s [ ' 
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to simplify the expressions. 

w. I ^ u - (B14) 

n^(M; 1 f) 1J - P art2 

^ ~ V (*+2)n M x - n U (k+2)n Am U (k+2)m ~ ~~ ^3 <-l (VLb) 

n^^O^-partS 

N W s ^ W ^ ^ *W - ^ / ° m, ' // ^' aili (B16) 
11^ (AC - Part 4 

^ ^ V (h+2)n U (P+2)n D ( k+2)m D (p+2)m ^ — ^ Cl (B17) 

n^^Oy-partS 



^ VS. I7 1B Am ^ 2)m - (B18) 



n^^Oy-parte 

^ ^ ft V2 " M - ^ Am 'W* - ft "l^f - Cl W s (B19) 

n^^O^-part? 

^ — V 2n M x - n U(h+2)n Am Ak+5)m - ~~ ^ C l ( B20 ) 

n^^Oy-partS 

-ft Jf 1 V (k+2)n M x - Uln V (k+5)m U {i+2)m ~ Jfi 2 <-l \P^) 



KpWn Oi,-part9 

* f 2 1 ^ M - ^ ^ * 5 >- * I ™ 2 C1 (B22) 

^^(M; 1 Ob - part 10 

v n/ - /X ' 7 ' A ' rj "' A ' 



w s 2 if v <*- 2 >» M - ^+ 5 >™ ^ +5 >- " "if mi 2 Ci (B23) 
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APPENDIX C: DETAILS OF (U^M^f) TERMS 

n° :- 

Note that the necessary symmetrization is not shown explicitly. 

• M Kkh ^'jhp ^d* _ o Kkh ^'jhp m d h {Ml L )* pk 

gT~gT dh VkmV ^ m ~ gTlh Ml ' 



s 



where 



2 V2M w cosB * V2M w smf3 
(-HJpfc = [A d - fx tan/3] m dp 6 kp + 5A kp ( 

( M 'n -P^t 1 

m ,, Ca ^ kh Vd h ™ r> d T> d * ~ -ihh Vdh Hi mdh {M^)hk or 



(AC -Part 2 

iu ,, ^4 A- fefc y dh , ,* A- fefc & h /ij m d „ {Ml L )* kh 

~ K N ~M s ~g~~g~2 dh hm ik+3)m ~~lhlh Ml ^ 



Ci XL 



Kkh to m *H™Q kh 3Ci 



tan^y M s 5> 2 



m d h T^km ^hm — 



92 



M% 



tan^ v 



Kb Parti 



Ci A' 



3M, 



£2 



m d h Aim Am — 



Kkh H m * h ™Q k 
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n Ss ^ Oy- Part 2 



~ 2 



2CiA- ftA , d d * A^ fe ^jin dh m 1 
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APPENDIX D: THE E„ RESULTS 



We have a simple result here, given as 
• 1 



[tan6W^ln - X 2n ] [£>(i +2 ) m - i £fi+7)m] [^0'+2)m + * ^(j+7)r 



BjB* tan 2 /3 



- (t&nOwXin — X 2n ) 2 



(n=l-4 dominate) . (Dl) 



We list all the individual terms below. 



^» D[ i+2)m D[ j+2)m ~ J (D2) 



~ 2 

- - U 2n U 2n D {l+b)m D (j+5)m _ - - ^ (D3) 
\ \* X A* 

— — ^+2)n ^(p+2)n U (k+S)m U (q+b)m ~ ~ — K^V 

i/2 i/2 i/2 i/2 

y " j tt tt* n r n l ~ iMi^hl (t)k\ 

— Uln U 2n D (i+2)m L> u+5)m ~ r-r^ (D5) 



9: 



2 



y 2rt - L/ (i+2)m ^(j+^m — ^ ^2 



■— U 2n U ln V {i+5)m V (j+2)m ~ — - — — (Db) 

i/2 i/2 JW S 



^hkfT TT* D l * D l ~ ^ fefc {^Ri)ki / n7 x 

— — L/ln ^(h+2)n ^(i+2)m ^(fc+5)m - ~ 7p l U 'J 

i/2 i/2 



^ihkjj jr* n l* n l ^ \hk Vh{M RL )%j / n8 x 

— — t/ ( h+2)„ C/ ln ^(fc+ 5 ) m ^(j + 2)m - — 7p \P%) 

i/2 i/2 JW S 

2n {h+2)n (t+5)m (fe+5)m - & 17 Ms (D9) 

~7 2 lh ih+2)n 2n (fc+5)m u+5)m -J^m's (D10) 
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V V V ^)n V U+ 2)n Am ^ ~ ~ *V ( D1 ^ 

i/2 i/2 i/2 i/2 

Aihfc , A i/tfc A* feg mg fcg + /4//q 

~ ~ V (A+2)n "(p+^n ^(fc+2)m ^( ? +2)m - ~ ~ TT 2 \- L>lZ > 

i/2 i/2 i/2 i/2 J|/J s 

i/2 i/2 i/2 i/2 

^ Aj feft j r Ajfcj g| tan/3 

" - — ^ + 2)n V (fc+2)B D (fe+2)m A m " ~ — (D14) 



92 92 92 92 

^Kkn^n * ^ ^ g Kkn ^jqn (m6) 

92 92 k 92 92 



APPENDIX E: SOME ILLUSTRATION ON THE VALIDITY OF THE 
APPROXIMATE FORMULAS THROUGH NUMERICAL CALCULATIONS 

In order to see how well our approximated formulas of the 1-loop neutrino mass correc- 
tions work, we present here some of the numerical neutrino mass values from the approx- 
imated formulas and compare them verses those from the exact expressions of the corre- 
sponding neutrino mass terms. We have a disclaimer to pronounce first. What we do here is 
not a numerical studies of phenomenological viable scenarios of neutrino masses generation 
within the model. We make no attempt to choose parameters to fit any neutrino oscillation 
data. Rather, we are choosing simple and quite arbitrary input parameters, only to check 
and give an idea on the validity of out analytical results. The practice also helps to illustrate 
some theoretical issues behind the formulas. We choose a set of convenient input parameters 
and compute and list results from nine of the long list of neutrino mass terms. While the 
results do give some idea on the relative strength of the various terms, the readers should 
be warned that this is only a consequence of a specific choice of inputs, which is in no sense 
generic or particularly phenomenologically interesting. 

Our choice of input parameters is as follows. We take the SUSY mass as around the 
scale of 100 GeV . In the exact results calculations, however, we have to split the masses 
of different superpartners to avoid unwanted special cancellations. We will clarify on the 
latter issue below. We choose input values that turned up mass eigenvalues for the SUSY 
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particles in the hundreds of GeV scale, details of which is really not interesting. The value 
of tan/3 is set at 3. The parameters responsible for the lepton number violating effect are 
simply taken to be the same numerically, at a value of 10~ 4 . Explicitly, = A' ijA , = 10~ 4 , 
Hi = 10~ 4 GeV, Bi = 10~ 4 GeV 2 . The neutrino mass results are presented in table 1, in 
which we show only contributions to the (3,3) elements of the effective (SM) neutrino mass 
matrix. The upper part of the table corresponds to the results from the approximated 
formulas while the lower part to those from the corresponding exact expressions. In each 
part, the first line corresponds to terns in Eqs.(42,43,B2), the second one to Eqs.(B3,B6,B15) 
and the third one to Eqs.(B22,48, CI). 

As one can see from the table, the difference between the two results for any specific 
contribution is within an order of magnitude. One could not really expect a better agreement 
than this. In fact, as discussed above and in some related earlier studies [17,18], the structure 
of the class of 1-loop diagrams are such that there is a GIM-type unitarity cancellation 
involved in the sum over mass eigenstates. Say, if all the mass eigenstate fermions of the 
same quantum number are degenerate, the sum over the set of fermion mass eigenstates in 
a neutrino mass term will be proportional to the mass matrix entry that a naive look at the 
Feynman diagram will suggest. In most cases, that is vanishing. Similarly, when the set 
of the mass eigenstate scalars involved in a certainly diagram is mass degenerate, the sum 
over the set of states gives a vanishing result due to unitarity of the diagonalizing matrix. 
Take Eq.(B6) as an illustrative explicit example, the exact expression of the neutrino mass 
contribution is proportional to 

In case of mass degeneracy, one can factor out a fermion summation 

E V (k+2)n M x - n U ^ 
n 

and a scalar summation 

LJ (k+2)m LJ (i+2)m ■ 

m 

The former is nothing but the vanishing (k + 2, 1) entry of the charged fermion mass matrix, 
while the latter is zero by unitarity of the matrix D l (for i ^ k). As discussed in Refs. [17,18], 
the lack of mass degeneracy leads to first order violation of such unitarity cancellations, which 
explains the nonvanishing results. It also explains the not better than order of magnitude 
agreement between our exact results, obtained really summing over all the contributions 
from the different mass eigenstates, and that from the approximate formulas, which only 
extract the analytical form of the largest term within such summations. 
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With the above explanation, we see that our approximate formulas do work as well as 
they are to be expected. We emphasize again that the approximate formulas mainly serve 
the purpose of illustrating the role of the lepton number violating parameters in each of the 
neutrino mass contribution term. 
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TABLES 

Table I. Some numerical results from the chosen neutrino mass terms. (See text of Appendix 
E). 

approximated formulae (eV) 



-3.85 x 1(T 9 -9.35 x 1(T 9 6.38 x 1(T 8 

1.78 x 1(T 10 -4.38 x 1(T 8 -2.47 x 1(T 9 

1.00 x 10~ 8 3.31 x 10~ 7 4.46 x 10~ 3 

exact formulae (eV) 

-2.61 x 10~ 9 -3.45 x 10~ 8 1.45 x 10~ 8 

5.63 x 10~ 10 -1.73 x 10~ 8 -3.94 x 10~ 9 

-3.23 x 10~ 9 6.90 x 10~ 7 3.26 x 10~ 3 
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